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ABSTRACT 

We show that the Hot-Dust-Poor (HDP) quasars, originally found in the X-ray 
selected XMM-COSMOS type 1 AGN sample, are just as common in two sam- 
ples selected at optical/infrared wavelengths: the Richards et al. Spitzer/SDSS 
sample (8.7% ± 2.2%), and the PG-quasar dominated sample of Elvis et al. 
(9.5% ± 5.0%). The properties of the HDP quasars in these two samples are 
consistent with the XMM-COSMOS sample, except that, at the 99% (~ 2.5cr) 
significance, a larger proportion of the HDP quasars in the Spitzer/SDSS sample 
have weak host galaxy contributions, probably due to the selection criteria used. 
Either the host-dust is destroyed (dynamically or by radiation), or is offset from 
the central black hole due to recoiling. Alternatively, the universality of HDP 
quasars in samples with different selection methods and the continuous distri- 
bution of dust covering factor in type 1 AGNs, suggest that the range of SEDs 
could be related to the range of tilts in warped fueling disks, as in the model of 
Lawrence and Elvis (2010), with HDP quasars having relatively small warps. 
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1. Introduction 

Emission from hot (~ 1500K) dust is so characteristic of AGNs and quasars (e.g. Sug- 
anuma et al. 2006) that the strong 1-3/zm emission from this dust has often been used to 
select AGN samples (e.g. Miley et al. 1985, Lacy et al. 2004, 2007, Stern et al. 2005, Donley 
et al. 2008). 

However, we recently reported (Hao et al. 2010) that, for 6% (at z < 2) to 20% (at 
2 < z < 3.5) of the quasars in the XMM-COSMOS type 1 AGN sample (Brusa et al. 2010; 
Elvis et al. 2011), their 1 — 3/im emission are two to four times weaker than typical type 
1 AGN (e.g., Elvis et al. 1994), although they have normal optical 'big-blue-bump' slopes 
indicative of standard accretion disk emission. We dubbed these 'hot-dust-poor' (HDP) 
quasars. Their implied dust 'torus' covering factor is ~2% to 30%, well below the 75% 
predicted by the unified model (e.g. Krolik & Begelman 1988). These sources lie at one 
extreme of the distribution of AGN SEDs found in the XMM-COSMOS sample. Their 
lack of hot dust emission uncovers a continuum that appears to be the continuation of the 
accretion disk, which implies a disk size of ~ 10 4 Schwarzschild radii, an order of magnitude 
beyond the gravitational instability radius. 

The study of Hao et al. (2010) used an X-ray selected sample. To check whether weak 
dust emission is a property of quasars in general, as opposed to X-ray selected quasars, in 
this paper we examine a large optical/infrared (IR) selected sample, that of Richards et 
al. (2006; hereafter R06). As shown in Figure [U R06 also covers a different part of the 
luminosity/redshift (L,z) plane from the COSMOS-XMM sample, which will help us to test 
the preliminary claim of Hao et al. (2010) that the HDP fraction changes with redshift but 
not with luminosity. We also examine the Elvis et al. (1994) sample of low-redshift quasars 
(hereafter E94). Figure [1] shows that the three samples together have a good coverage of the 
(L,z) plane. 

2. Type 1 AGN samples 

Hao et al. (2010) selected HDP quasars by their position in a plot of optical (OPT, 0.3- 
1/xm) versus near-infrared (NIR, 1-3/im) slopes. In this paper, we use the samples of R06 
and E94 supplemented by new near-IR UKIDSS data, which gives good enough photometry 
to define these slopes. 
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Fig. 1. — The bolometric luminosity and redshift of the XMM-COSMOS (squares), R06 
(circles) and E94 (triangles) sample. The red solid points are the HDP quasars from these 
samples. The E94 quasars show a narrow distribution, as they are chosen to be the very 
brightest quasars in the sky. is calculated as in Hao et al. (2010) by integrating from 
rest frame 24 [im to the Lyman limit (9 12 A) and using the WMAP 5- year cosmology: 
H = 71 km s^Mpc' 1 , Q M = 0.26 and Q A = 0.74 (Komatsu et al, 2009). 
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2.1. R06 sample 

The R06 sample consists of 259 Spitzer sources identified with SDSS quasars in four 
discrete degree-scale fields, and is, therefore, mid-IR and optically selected. We also consid- 
ered the similar Hatziminaoglou et al. (2008) sample of 278 quasars and found that all but 
65 are also in R06 and only 39 of these 65 have Spitzer photometry. We therefore only use 
R06. The redshift range covered is z = 0.14 — 5.2 with 93% being at z < 3 (Figured]). 

Most of the data in R06 paper is from SDSS, Spitzer, and 2MASS photometry. These 
data were not corrected for host galaxy contributions. As R06 is an IR and optically selected 
sample, the host galaxy contribution should typically be small. As the galaxy SED in the 
IR and optical is completely different from the typical quasar SED, and IR and optical color 
selection picks out sources with spectral indices similar to those of a typical quasar, sources 
with more galaxy contribution are selected against (e.g. Gregg et al. 2002, Richards et al. 
2003). Most (215/259) of the R06 sources did not have 2 MASS J H K photometry. These 
bands are essential to make reliable NIR slope estimates for quasars at z < 0.65, and to 
make reliable OPT slope estimates for quasars at z > 1.3. We therefore cross-matched the 
R06 sample with the UKIDSS database using the online WFCAM science archive. We found 
additional data for 98 quasars at K and for 45 quasars at J, producing a final subset of 195 
R06 quasars that we can use for this study. 

The UKIDSS data come from the Deep Extragalactic Survey (DXS) fields covering the 
Lockman Hole and ELAIS-N1. The UKIDSS survey is defined in Lawrence et al (2007). It 
uses the UKIRT Wide Field Camera (WFCAM; Casali et al, 2007). The photometric system 
is described in Hewett et al (2006), and the calibration is described in Hodgkin et al. (2009). 
The pipeline processing and science archive are described in Irwin et al (2010, in prep) and 
Hambly et al (2008). The magnitudes we use here are "aper3" magnitudes, which come 
from a measurement through a 2 " software aperture, corrected to a total magnitude using 
the point spread function measured in the field. We have not corrected these magnitudes 
for host galaxy contribution, but we could see that nearly all of these objects only have 
small extended source contribution, because the UKIDSS database provides a wide range 
of aperture magnitudes. The magnitudes we use come from stacked observations spreading 
over several years. From individual epochs, we can see that most of these quasars vary only 
by a small amount (0.05 mag) within the UKIDSS data, such that it does not seriously affect 
our analysis. 
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2.2. E94 sample 

The E94 quasar sample consists of both optically-selected quasars from the bright quasar 
survey (BQS, PG, Schmidt & Green, 1983) and radio-selected (mainly 3C and PKS) quasars. 
These quasars were selected to have good signal-to-noise ratio Einstein X-ray observations 
and obtainable IUE UV spectra, hence E94 is biased towards X-ray bright and blue quasars. 
E94 contains 42 quasars in the redshift range z = 0.025 — 0.94, with 80% being at z < 
0.3. The optical photometry were obtained at the FLWO (F. L. Whipple Observatory) 
24 inch telescope within one week of the MMT FOGS (Faint Object Grism Spectrograph) 
observations. The NIR data were obtained with MMT and IRTF. More details on the 
observation can be found in Elvis et al. (1994). As E94 is mainly a local quasar sample, 
the E94 SEDs have been corrected for host galaxy contamination by subtracting the host 
galaxy template SED based on the Sbc galaxy model of Coleman, Wu & Weedman (1980). 
The normalization of the galaxy template was calculated by direct measurements of the host 
galaxy luminosities. 



3. Mixing Diagram 

To select HDP quasars, we use the mixing diagram technique of Hao et al. (2010), 
where we select objects that lie outside the region formed by mixing a typical quasar SED 
with a host galaxy SED, combined with reddening. We use the slope on either side of the 
1/zm inflection point, which is where the black-body emission of the hottest dust rises above 
that of the accretion disk. We build the mixing diagram by plotting two rest frame slopes in 
log vL v versus log v space: ctNm (3/ira — lfim) versus cxopt (lfim — 3000A). These slopes 
are derived from linear fits to the SEDs in these wavelength ranges. The errors on the slopes 
are the la standard error. We required that at least three photometry points are used to fit 
the slopes. 88% of R06 quasars have more than the minimum of three photometry points 
in the optical and 71% have more than three points in the near-infrared. The measurement 
error on the photometry is used in the fitting. For both R06 and E94, all these slopes are 
measured in the same rest frame wavelength range as the XMM-COSMOS sample. 

The mixing diagram is equivalent to a color-color plot, but utilizes more photometric 
points and adapts well to a wide range of redshifts. The mixing diagram is more robust 
than the normal color-color plot, because it is less affected by one bad photometry point and 
is built in the rest frame without any need for assumptions about the intrinsic SED shape 
(k-correction) . 



-6 - 



3.1. Robustness of Linear Fitting 

We performed two tests to check if the requirement to have three or more photometry 
band is adequate to ensure robust slope measurements. 

(1) We use the E94 mean SED as a template. This SED has aopr = 0.95 and a^jp = 
—0.72. We made a Monte Carlo sampling of the E94 SED, with the same sampling and 
error bars as the observed SEDs for each R06 quasar. We then measured the simulated 
slopes, as we did for the observed ones. We show the histogram of the simulated slopes 
in Figure EJ 

For the simulated optical slopes, the mean is {cxopt,mc,E9a) = 1-07 (<Jopt,mc,E94 = 
0.19). The E94 optical slope (ccopt = 0.95) is within 1ctopt,mc,e94 of the mean value. 
However, when considering the errors (a) on each Monte Carlo optical slope, we find 
that only 75% of the simulated quasars lie within 3a of the E94 optical slope, implying 
significant outliers. In fact, a tail toward steep slopes is evident in Figure [2] (left). 
This tail is cause by the curvature in the E94 SED in the 0.3-1/zm range. The steepest 
part of E94 reaches a slope of ~ 1.5 between 0.3-0.45/iin (14.8 < logu < 15). If the 
simulated photometry points happen to lie only in this region, then the simulated slope 
will have a larger value than the o>opt = 0.95 calculated over the full 0.3-1/zm range. 
This is a limitation of using just the E94 template, which we will address next. 

For the corresponding simulated infrared slopes, the mean is ((Xnir,mc,es±I = —0.70 
(&nir,mc,eh = 0.17). The E94 infrared slope (a NIR = -0.72) is within 1(t N ir,mc,E9A 
of the mean value. When considering the errors (a) on each Monte Carlo infrared slope, 
we find that 91% of the simulated quasars lie within 3a of the E94 infrared slope. The 
histogram of the simulated infrared slopes shows two peaks. One peak is located at 
~ —0.72, which is the E94 infrared slope. The other peak is located at ~ —0.85. As 
for the tail in cxopt,mc,E9^ this second peak is due to curvature in the E94 SED in 
the 1-3/zm range. The steepest part of E94 reaches a slope of ~ —0.9 at 1.2-2.4/xm 
(14.1 < logv < 14.4). If the simulated photometry points happen to lie only in this 
region, the simulated slope will have a large value compared to the E94 infrared slope 
calculated over the full 1-3/im range. 

While the mean of the simulated slopes are same as the E94 values, the offsets are a cause 
of concern as they may falsely indicate an HDP object. Therefore, we performed a second 
Monte Carlo test using the measured slopes of the COSMOS sources with good photometry 
coverage. This method avoids the systematic errors from using a single template SED. 
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Fig. 2. — Histogram of the E94 slopes with the SDSS-Spitzer sampling for the R06 quasars: 
Left: oiopr] Right: a^iR- The blue dashed lines are the Gaussian fit of the histogram. 




Fig. 3. — Comparison of the simulated SDSS photometry slopes with the XMM-COSMOS 
photometry slopes for the XMM-COSMOS quasars: Left: a opt', Right: onir- The red solid 
line is the equal slope line. The green line is the linear fit of the black points. The two 
red dashed lines are at the 3(7^ distance to the equal slope line. The error bars of the slope 
measurement are shown as yellow. 
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(2) To test if the slope measurement is robust for the range of SED shapes seen in the XMM- 
COSMOS type 1 AGN sample, we used the XMM-COSMOS SEDs as the template 
SEDs. The photometry coverage for these SEDs is rich in the NIR to OPT range, 
giving robust slope measurements (Elvis et al. 2011, Hao et al. 2010). 

We selected 398 (97%) XMM-COSMOS type 1 AGNs, that lie within Az = 0.1 of a 
corresponding quasar in R06 sample. We used these 398 quasars to perform a Monte 
Carlo sampling, matching the photometry and error bars to those observed in the 
corresponding R06 quasar. We then measured the resulting ocopt,mc and aNm,MC, 
and compared them with the original measurements from COSMOS. The results are 
shown in Figure [31 A linear fit to aopT,MC versus «opt nas a slope of 0.98±0.0004. 
The linear fit of &nir,mc versus a^m has a slope of 1.005±0.0003. Both slopes are 
very close to 1. However there are outliers that emulate HDP quasars. 

To investigate these outliers, we defined the distance of each data point to the equal 
slope line as d w = (a W; MC ~ «w)/v / 2 where w = OPT, NIR. A positive value means 
that the point lies above the equal slope line. For the simulated optical slopes, the 
mean of dopr is 0.015 (u^opt = 0.20), consistent with in Ict^opt- For the simulated 
near-infrared slopes, the mean of d^m is 0.067 (cr d>OPT = 0.15), consistent with in 

lCd,7V_TR- 

We use these cr^w to isolate 16 sources that lie outside the "io^NiR lines (red dashed 
lines), 4% of the total sample. These sources have 3 points in the NIR slope estimation 
and at redshift 0.92 < z < 1.3 (see example in Figure El right). We plotted the d^ip 
versus redshift, shown in Figure HI For an E94-like SED shape, the three IRAC bands 
in this redshift range would lie at > 2/im, in the upper part of the infrared bump. The 
lack of a K band measurement would then lead to a larger slope estimate. For galaxy- 
like SED shape in this redshift range, lack of K band would lead to a much smaller 
slope estimate. For a mixing of E94 and galaxy SED, which is quite flat, the lacking of 
K band will not affect the slope measurement. The total number of simulated sources 
with 3 points in the NIR slope estimation is 95. Thus for the majority (83%) of the 
sources with 3 photometry points in the NIR slope estimation, the slope estimates are 
still robust. 

In all, for the 96% simulated sources with 3 or more photometry points in NIR, the 
slope estimation is robust. The NIR slope estimate is not quite reliable for sources 
with only 3 photometry points and at redshift 0.92 < z < 1.3. 

In general, for majority of quasars in R06 sample, the slope measurement we applied 
is robust for E94 like SED shape or a variety of SED shapes as in XMM-COSMOS sample, 
even for quasars with limited photometry coverage as in R06. Figure H] also shows that 
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the systematic error on at^m is 0.07 for sources using 4 or more photometry points, which 
is small compared to the measurment errors. For sources using only 3 photometry points 
in the redshift range 0.92-1.3, the systematic error on olnir is 0.3, which is comparable in 
some cases to the measurement errors. We added this errors to the measurement error when 
plotting the mixing diagram of R06. 

3.2. Corrections to the SED 

We need to consider the following several factors, which might affect the measurement 
of the slopes. 

3.2.1. Emission Lines 

The major emission lines in the rest frame wavelength range 0.3 — 3/im are the Balmer 
series, the Paschen series, and the Fe II complex. As the photometry points in the SED are 
all broad-band, even for the strongest line, Ha (mean equivalent width 257A for SDSS DR6 
quasars, Elvis et al. 2011), the upward bias of the observed photometry from the continuum 
reaches only 0.1 dex in the z band (filter FWHM=1000A). For weaker lines, such as H/3 and 
Fe II, the bias would be ~0.02 dex, which is negligible for the slope measurements reported 
here. For consistency, we correct the R06 photometry for the emission lines as we did for 
the XMM-COSMOS sample (Elvis et al. 2011). 

For the E94 sample, the emission line contribution has already been corrected for in the 
reported SEDs (details in Elvis et al. 1994). 

3.2.2. Variability 

Most AGNs vary in their optical continuum flux by ~10% on timescales of months to 
years (Vanden Berk et al. 2004). The photometry used to build the SEDs in the R06 sample 
were taken over different time periods, which may affect the estimation of slopes. 

For the R06 sample, the IRAC detections were made around 2004 (Lacy et al. 2005, 
Surace et al. 2005). The SDSS-DR3 photometry including the optical SDSS bands were taken 
from 1998 to 2003 and the NIR 2MASS bands were taken from 1997 to 2001 (Schneider et 
al. 2005). The UKIDSS NIR photometry has been taken since 2005. We should be careful 
in using the UKIDSS data and the rest of the data to measure the slopes. 
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For the olnir estimation, we need to consider sources in the following redshift ranges 
respectively: 

1. z > 1.15: only the IRAC bands lie in the 1-3/im range. Thus the a^m estimation is 
robust, for the 109 out of the 195 R06 quasars; 

2. 0.65 < z < 1.15: K and IRAC bands lie in the l-3//m range. 16 sources have UKIDSS 
K-band detection. In no case does the K-band point appear significantly discontinuous 
with the IRAC bands. 

3. 0.24 < z < 0.65: H, K and IRAC bands lie in the 1-3/im range. 17 sources have 
2MASS H band detections. 8 sources have UKIDSS K-band detection. In no case does 
the H-band and K-band point appear significant discontinuous with the IRAC bands. 

4. z < 0.24: All of the J, H, K and IRAC1 bands lie in the 1-3/mi range. 2 of the 3 
sources at z < 0.24 have UKIDSS J and K band detections. The SEDs of these 2 
sources show no discontinuity in the SED to cause errors in the slope measurement. 

We thus consider the estimates of a^m to be robust. 

For aopT estimation, similarly, the 3 low redshift (z < 0.24) sources have only the 
SDSS data in the 0.3-1/iin range. Thus for these quasars the «opt estimation is robust. 
For sources at 0.24 < z < 1.15, only the J band and SDSS data lie in the range. 12 sources 
have UKIDSS J band detection, for which the effect of J band on the slope measurement 
is small. For the 109 high redshift (z > 1.15) sources, the situation is more complicated, 
because the UKIDSS J K band and the SDSS or 2MASS H band could be in the 0.3-1/zm 
range. However, the effect of variability on the aopr does not affect the selection of the HDP 
quasars, since the most significant difference of the HDP quasars from the normal quasars 
are the a^m- 

For the E94 sample, the optical and the ground-based IR data were generally obtained 
within about 1 month (Elvis et al. 1994), so that the variability does not affect the slope 
measurement in the optical and IR range used in the mixing diagram. 

3.3. R06 and E94 Mixing Diagram 

Figure [5] shows the Hao et al. (2010) "mixing diagram" for the R06 and E94 samples. 
For the sources that have galaxy fraction larger than 0.6 in the mixing diagram, we check 
the inferred host galaxy luminosity. Only one quasar in R06 sample has M < —23 (-23.8) 
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at 1/xm. Thus the inferred host luminosity is reasonable. There is clearly a continuous 
distribution of objects in a^jn, so to select quasars in the HDP region, we need to examine 
the dispersion in "normal" objects. 

We consider all R06 quasars with 0.2 < aopT < 1-6 in order to exclude the few galaxy- 
or reddening-dominated sources. They have mean slopes olopt = 0.76 (standard deviation 
Copt = 0.25), and cunir = —0.53 (standard deviation ct^ir = 0.37). The mean of the 
measurement error of olopt is Ettqpt — 0.10, which indicates the intrinsic dispersion is 
&int,opt — a/ o"opt — Err OPT = 0.23. The mean of the measurement error of olnir is 
Err^jji = 0.08, which indicates the intrinsic dispersion is &int,nir — a/ crf fIR — Err 2 NIR = 
0.36. The intrinsic dispersion of the SED shape is substantial. 

The E94 sources, which were corrected for host galaxy contribution, are, by construc- 
tion, clustered around the E94 mean in the AGN-dominated region located at the bot- 
tom right corner of the mixing diagram. For the 42 E94 sources, the mean slopes are 
otoPT = 0.96 (standard deviation oopt = 0.27) and cxnir = —0.50 (standard deviation 
&nir = 0.38). The mean of the measurement error of aopr is Err opt = 0.1, which in- 
dicates the intrinsic dispersion is &int,opt = a/ o~ opt — Err OPT = 0.25. The mean of the 
measurement error of a^m is Err NIR = 0.2, which indicates the intrinsic dispersion is 

VINT,NIR = V a NIR- Err NIR = °- 32 - 

To compare the intrinsic dispersion of these two samples, we performed an F-test, ig- 
noring the effects of the measurement error and assuming the intrinsic slope distribution 
is normally distributed. For aopT, the F statistic is 1.09, and the cumulative distribution 
function (CDF) probablity is 65%, implying that the difference between the olopt intrinsic 
dispersions of R06 and E94 is not significant. For cxnir, the F statistic is 1.13, the CDF 
probability is 66%, implying the difference between the qinir intrinsic dispersions of R06 
and E94 is not significant. A more detailed statistical comparison between the intrinsic 
dispersion in different samples is reported in Hao et al. (2011). 



3.4. Selection of the HDP quasars 

In the mixing diagram, we define a circle with a radius of 0.6 to approximate the 1. So- 
intrinsic dispersion region of the E94 mean SED template. We then use this region to classify 
outliers. As the distribution of the quasars is continuous in a^m, different sized circles define 
different outlier populations. Here we simply apply the same criterion as in Hao et al. (2010) 
to enable ready comparison of their properties. 

The Hao et al. (2010) mixing diagram readily distinguishes among the AGN-dominated, 
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galaxy-dominated and reddening-dominated SEDs. The AGNs with SED slopes in the tri- 
angular region defined by the mixing curve and the reddening vector can be explained by 
the combinations of the E94 mean SED, host galaxy contamination and reddening. AGNs 
which lie at least la above the mixing curve and the dispersion circle in the upper right 
corner of the mixing diagram, are defined as HDP quasars. Note that we do not suggest 
that HDP quasars are a distinct class. As the distribution of the quasars on the mixing 
diagram is continuous and the definition of the dispersion circle is somewhat arbitrary, the 
HDP quasars represent one extreme of the distribution of the AGN SEDs. 

We further divide the HDP quasars into three classes by their positions relative to the 
equal slope line as in Hao et al. (2010). Class I sources (blue points) have aopr > &nir 
within measurement error. These sources have a normal blue bump but flat NIR emission. 
Class II sources (green points) lie on the equal slope line, which indicates that the NIR 
emission might be the continuation of the optical accretion disk, although this turns out to 
be not true for all cases judging from SED fitting in the XMM-COSMOS sample (Hao et al. 
2010). 

We find 17 HDP quasars in the R06 sample. Most (15) of which belong to Class I, 
with the remainder (2) belonging to Class II. The fraction of the HDP quasars in R06 is 
8.7% ± 2.2%. For 3/4 of these HDP quasars, more than 3 photometry points are used in 
the near-infrared for slope estimation. The two Class II HDP quasars, however, have only 3 
NIR data points in their estimates. 

Four sources out of the 42 quasars in the E94 sample qualify as HDP (Table [2]): 2 
optically selected and 2 radio selected. Their SEDs are shown in Figure [HJ They are all 
Class I HDP quasars. 

4. HDP quasar properties 

4.1. HDP Quasars in R06 

Representative R06 SEDs (Class I: top, Class II: bottom) are shown in Figure |6J The 
Class I source (SDSS J104114.18+590219.4) by definition show a flat near infrared SED, 
similar to the Class I HDP quasars found in XMM-COSMOS quasar sample (Hao et al. 
2010). 

The Class II source (SDSS J104935. 76+554950.5) has only 3 points in both slope esti- 
mates and is at z = 1.056, in the range where a^jR can be overestimated by up to ~ 0.5 
(Figure H] in section |3TT|) . However it is hard to fit the SED with an E94 mean SED and some 
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Fig. 4. — The distance of the near-infrared slopes (Figure [3j right) to the equal slope line 
versus redshift. The red squares represent the Monte Carlo slope measured with 3 points. 
The black cross represent the Monte Carlo slope measured with more than 3 points 

Table 1. HDP Quasars in R06 



Name (SDSS J) redshift logL boJ Class Afjvii? a 

erg-s — 1 



105308.24+590522.2 


0.430 


45.16 I 


4 


160341.44+541501.5 


0.581 


45.27 I 


4 


104857.92+560112.3 


0.800 


45.70 I 


3 


103525.05+580335.6 


0.964 


46.02 1 


4 


163306.12+401747.4 


0.974 


45.87 I 


4 


160913.18+535429.6 


0.992 


46.32 1 


4 


104935.76+554950.5 


1.056 


46.26 11 


3 b 


104114.18+590219.4 


1.093 


45.93 1 


4 


104921.49+575036.6 


1.106 


46.08 1 


4 


171441.11+601342.9 


1.456 


46.05 1 


4 


104948.86+592620.7 


2.000 


46.22 I 


4 


103905.24+591209.1 


2.115 


46.50 I 


4 


104226.20+585925.8 


2.560 


46.48 I 


4 


160828.33+535251.9 


2.960 


46.61 I] 


3 


164238.08+412104.7 


3.123 


46.76 1 


4 


104810.81+575526.9 


3.171 


46.71 1 


4 


105654.96+583712.4 


3.856 


46.81 I 


3 



a Number of photometry points used in the «jv/fl estimates. 

h &NlR may be overestimated (see section [3.11 Figure |4). Detailed discussion about this source is in section [4+] 
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Fig. 5. — Mixing Diagram: Left: R06 sample; Right: E94 sample. Red cross and red circle 
show the E94 mean SED and the la dispersion of the E94 sources. The blue triangles show 
16 different SWIRE galaxy templates (Polletta et al. 2007). The purple lines connecting the 
E94 and the galaxy templates are mixing curves showing the slopes of different fraction of 
galaxy and AGN. The straight purple arrow shows the reddening vector of E94. The straight 
red solid line shows the oiopt — oinir line. Different colors of the points show different class 
of the HDP sources (I-blue, II-green, see text for details). The black symbols show all 
the other type 1 AGN in the samples. For R06, the error bars are adjusted for possible 
systematic errors due to limit photometry. The 6 magenta dots in the mixing diagram of 
R06 sample show the z ~ 6 quasars in Jiang et al. (2010) Figure 1. 
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host. The NIR SED of this source shows an unusual discontinuous, 0.1 dex arise from the 
IRAC2 to IRAC1 points at ~ 2/xm (rest). This feature cannot be caused by the emission 
from hot dust, because a black-body, at any temperature below the maximum dust subli- 
mation temperature of 1900 K, does not fit the feature well. The Paschen a line lies in the 
IRAC 1 band but is unlikely to produce this feature, as to produce the IRAC1 rise would 
require an equivalent width of ~2000A, which is 30 times the typical value of 70A (Landt et 
al. 2008). As «opt is within the range of those for E94 quasars (section I4~2"|) . it could be that 
the IRAC1 point is dominated by the continuation of the accretion disk. NIR photometry 
or spectra are needed to understand the emission in this range. 

As reported in Hao et al. (2010), the fraction of the HDP quasars in the XMM-COSMOS 
sample evolves with redshift but not with other parameters (e.g., Lw). The fraction increases 
from 6% at low redshift (z < 2) to 20% at moderate high redshift (2 < z < 3.5). We checked 
if similar evolution behavior can be found in the R06 sample. We plotted the HDP fraction 
in R06 versus z and Lboi in Figure [71 The HDP quasars in the R06 sample have redshifts from 
0.43 to 3.86. The range of L bol for the HDP AGNs (45.1<logL feo , <46.8) is almost the same 
as for the whole sample (44.8<logLfe z <47.4). The Kolmogorov-Smirnov (KS) test shows 
that the distributions of z and for the R06 HDP quasars are indistinguishable from the 
rest of the R06 sample (KS probabilities of 0.25 and 0.88). As the XMM-COSMOS sample, 
the HDP quasars in R06 do not show a different Lboi distribution to the other quasars in R06. 
However, the HDP quasars in R06 do not show redshift evolution, unlike XMM-COSMOS 
sample. This result might be affected by the 64 sources we have excluded due to limited 
photometry coverage. At redshifts z < 0.5, the 5 sources with no J H K photometry have 
fewer than three photometry points in the rest frame 1/xm — 3/im band, as do the 6 quasars 
at high redshifts z > 3.7. At redshifts 1.5 < z ^ 3.7, 53 quasars also have fewer than 
three photometry points in the rest frame 0.3 — \[im. So quasars at low and high redshifts 
have more chance of being excluded from consideration, which might lead to the lack of 
an increase in HDP fraction in R06 at z ~ 2, though «opt only has a small influence on 
HDP definition. The XMM-COSMOS sample has no such incompleteness bias, since the 
COSMOS photometry coverage is more complete and uniform. 

In the XMM-COSMOS sample, the HDP fraction jumps at z = 2 from ~ 6% to ~ 20%. 
So we re-binned the R06 sample and find that the HDP fraction is 6.3% ±2.1% at z < 2, and 
19.4% ± 8.0% at z ^ 2 respectively, a 1.6<r difference, which is consistent with the results in 
the XMM-COSMOS sample. We over-plot XMM-COSMOS fit of the HDP fraction versus z 
(red lines, f HDP = (0.031 ± 0.002)(1 + z)( L2±01 ), Hao et al. 2010). The HDP fraction versus 
z relationship in the two samples do not disagree. 

Because of the poorly sampled SED at ~ 1/xm due to limited NIR photometry in the 
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R06 sample, we could not make robust estimates of the dust covering factor (f c ) or the outer 
edge of the accretion disk (r out ). 

4.2. HDP Quasars in E94 

The fraction of HDP AGNs in E94 is 9.5% ± 5.0%. Because the sample size of E94 is 
10 times smaller than that of the XMM-COSMOS type 1 AGN sample, the error on the 
fraction is much larger. A KS test shows that the distributions of z and for these HDP 
quasars are indistinguishable from the rest of the E94 sample (with KS probabilities of 0.5 
and 0.6 respectively). These results agrees with the XMM-COSMOS sample that for low 
redshift quasars, the fraction of the HDP do not evolve with z and Lw 

As the E94 SEDs are all corrected for host galaxy contributions, we fit the SED with 
just two components: an accretion disk (using a standard Schwarzschild a-disk model, with 
electron scattering and Comptonization; Siemiginowska et al. 1995), and a hot dust compo- 
nent (using a single temperature black-body). The black hole masses for these four quasars, 
reported in Vestergaard & Peterson (2006), Marziani et al. (2010) and Woo et al. (2002), 
is listed in Table |2j Their accretion rates (expressed as their Eddington ratios, Xe) are cal- 
culated from Lboi and Mbh- We use the accretion disk model (Siemiginowska et al. 1995) 
for the particular black hole masses and the accretion rates for each quasar and show the 
results in Figure [HJ 

4-2.1. Covering Factors of E94 HDP Quasars 

Following the same method as in Hao et al. (2010), the dust covering factor (f c ) was 
estimated as the ratio of the dust emission area (Ad) to the area at the dust evaporation 
radius (A e = 47rrg). Ad comes from the normalization of the black-body fit. The dust 
evaporation radius (r e ) is estimated from r e = I.SLI/^qT^IqPc (Barvainis 1987), where 
L U v,46 is the total ultraviolet (l/xm-912A) luminosity in units of 10 46 erg s -1 , and T 150 o is the 
maximum dust temperature allowed by the SED, in units of 1500 K. The typical errors on 
the E94 photometry in the near- infrared are ~ 0.1 dex. This produces an error on the dust 
temperature estimation, ATd/Td ~ 8%, and an error on the normalization, A Ad /Ad ~ 4%. 
The error on f c is thus Af c / f c ~ 5%. 
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Fig. 6. — Examples of HDP SEDs in R06: Left: Class I; Right: Class II (aopr and olnir 
less reliable see section I37T]) . The red dashed line is the E94 radio quiet mean SED. 



Table 2. HDP Quasars in E94 



Object 


Name 


rcdshift 




log(M sff ) 




T d 


A d 


fc 


Tout 


Rout 


Rout/ 


Rout/ 








erg-s - 1 


M Q 




K 


pc 2 


% 


K 


pc 




T grav 


Q0003+158 


PHL 658 


0.45 


46.5 


9.3 1 


0.13 


1200 


8.00 


6 


3700 


0.6 


3200 


11 


Q0049+171 


PG 0049+171 


0.064 


44.8 


8.3 1 


0.03 


1300 


0.13 


11 


4000 


0.06 


2700 


11 


Q0414-060 


3C110 


0.78 


47.0 


9.9 2 


0.10 


1900 


3.80 


18 


3600 


0.4 


370 


13 


Q1635+119 


MC2 1635+119 


0.146 


45.4 


8.1 3 


0.16 


1900 


0.15 


21 


4800 


0.08 


6300 


9 



1 Reference: Vestergaard & Peterson (2006) 
2 Reference: Marziani et al. (2010) 
3 Reference: Woo et al. (2002) 
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Fig. 7. — The redshift distribution (top) and fraction of the HDP AGN in R06 sample 
as function of redshift (middle) and (bottom). In the top panel, the black histogram 
shows the redshift distribution of the R06 sample. The blue shaded histogram shows the 
distribution of the HDP quasar in R06 sample. The two histogram are scaled to have the 
same integral area. In the middle panel, the red solid line shows the fitting of the XMM- 
COSMOS HDP fraction versus redshift in the same bin (f HDP = (0.031±0.002)(l+z)( L2±a1 )) 
and the red dashed line show the la range of the fitting. 
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Fig. 8.— HDP AGN in E94 sample: Left: Q0003+158 (PHL 658, redshift z=0.45); Middle 
left: Q0049+171 (PG0049+171, z=0.064); Middle right: Q0414-060(3C110, z=0.78); Right: 
Q1635+119(MC2 1635+119, z=0.146). The SEDs are fitted with the accretion disk compo- 
nent (red dashed line) and a hot dust component (magenta dashed line). The sum of the 
two component is showed as black dotted line. The blue dashed line (a single temperature 
black body) shows the fitting of the outer edge of the accretion disk. The green dot-dashed 
line is the spectrum from a steady optically thick accretion disk (uF u oc z/ 4//3 ). 
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4-2.2. Accretion Disk Outer Radii of E94 HDP Quasars 

The outer accretion disk radius is estimated from the standard a-disk model formula 
(Frank, King & Raine, 2002) 

R out = 1.1 x 10 4 TPa~^r]--sMf 5 \ifl. pc 



where M 8 = M/(10 8 M o ), X E = L bol /{^^M), and / = [l - (ff ) 
a = 0.1, r] = 0.1. Thus R out , pc (x T ~f M^A| ex T~f Ms. 



We assume 



The error in the disk temperature estimates from the E94 photometric errors is /\T ou t /T ou t 
5%. Unfortunately, black hole mass estimates from mass scaling relationships have an error 
AM BH /M BH ~ 40% (Vestergaard & Peterson 2006, Peterson 2010), so AR out<pc /R outtPC ~ 
20%. In terms of Schwarzschild radii (r s = 2GM/c 2 ), AR out ,r a / Rout,r s ~ 33%. 

We can also express this radius in terms of the gravitational instability radius of the 
accretion disk (r grav , Goodman 2003), which is the radius beyond which the disk is unstable 

_2 ^ _2 -4 

to self-gravity and should break up. As r grav oc M sAJ, oc M 3, R ou t, rgrav oc T ou l . Hence, 

ARout,r grav / Rout,r srav ~ 7%. 

The results are reported in Table El They turn out to be comparable to the HDP quasars 
in the XMM-COSMOS sample in that R out ~ 10r grav , a surprisingly large value. 



5. Discussion and Conclusions 

5.1. Comparison with COSMOS HDP Quasars 

We find a similar fraction of HDP quasars in the optically selected R06 and E94 samples 
as in the X-ray selected XMM-COSMOS sample of Hao et al. (2010): 8.7% ± 2.2% in the 
R06 sample and 9.5% ± 5.0% in E94 sample. This similarity indicates that these extreme 
AGNs are fairly common even in optically selected samples, though they were not previously 
recognized. 

The XMM-COSMOS quasar sample and the R06 sample have similar redshift coverage. 
In contrast to the XMM-COSMOS sample, however, the fraction of the HDP in R06 does 
not show significant evolution with redshift according to the KS test result. This difference 
might be caused by the exclusion of the 64 sources in certain redshift range due to the limit 
near-IR photometry. The evolution of the HDP fraction in R06 sample is ill-determined but 
generally agrees with the HDP of XMM-COSMOS sample (Figure (H). 
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In the small E94 sample at low z, the redshifts and luminosities of the HDP quasars are 
indistinguishable from the normal type 1 AGNs. This is consistent with the XMM-COSMOS 
sample, as for low redshift sources {z < 2), the HDP fraction does not evolve with redshift 
and luminosity. 

All of the HDP quasars in R06 are Class I and II. In comparison, the XMM-COSMOS 
sample has 27% (11 out of 41) Class III HDP quasars. We compare the Class III distribution 
versus Class I and II in the R06 HDP and XMM-COSMOS HDP samples using the Fisher 
exact test (Wall & Jenkins 2003). The probability that the two HDP samples have the 
same class distribution is 1.4%. The SEDs of the Class III HDP quasars have stronger 
galaxy contribution than Class I and II. So R06 includes fewer sources with strong galaxy 
contribution than XMM-COSMOS. This is expected, as the R06 sample is selected by optical 
colors, and will exclude low AGN-to-galaxy contrast objects. The XMM-COSMOS sample, 
instead, is X-ray selected, and includes sources with strong galaxy contributions. 

In E94, the dust covering factor of the HDP quasars ranges from 6% to 21%, similar 
to the XMM-COSMOS HDP quasars (Hao et al. 2010). A small dust covering factor 
was proposed long ago in several quasars. For example, the spectropolarimetry analysis of 
the quasar OI 287 showed that the broad Balmer lines are polarized the same way as the 
continuum, while the forbidden lines have zero intrinsic polarization, suggesting a thin 'torus' 
(Goodrich & Miller, 1988). 

The outer edge of the accretion disk in E94 (Table [2]) ranges from 0.06pc to 0.6pc and 
from 370r s to 6300r s , which are ~ 10 times greater than the gravitational instability radius 
(Goodman 2003). A similar extension of the disk spectrum into the NIR uncovered by 
polarized light, was found by Kishimoto et al. (2008). The outer radius of the accretion disk 
is further out than expected by simple a-disk theory. For the R06 sample, the limitations of 
the photometry prevents the accurate estimation of these quantities. 

5.2. HDP Quasars and IR Quasar Selection 

We checked whether the HDP quasars lie within the Lacy et al. (2004) color-selection 
region for quasars (Figure [9]). Most of the HDP quasars lie in the bottom left of the quasar 
region. Although their hot dust emission is weak, it is still sufficient to differentiate them 
from starburst galaxies. The HDP quasars are one of the extreme of the variety of the AGN 
SEDs. AGNs with galaxy-dominated SEDs lie in the same region. The Lacy et al. (2004) 
color-color plot does not distinguish the HDP quasars from galaxy-dominated quasars. 
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5.3. Comparison with Jiang et al. Hot-Dust-Free Quasars 

Jiang et al. (2010) introduced the rest-frame 3.5/jm to 5100A luminosity ratio as an 
indicator of the hot-dust abundance. They concluded that there are no low redshift counter- 
parts to the two z ~ 6 hot-dust-free quasars J0005-0006 and J0303-0019. We reproduce the 
Jiang et al. (2010) plot in Figure dS for the XMM-COSMOS, R06 and E94 samples. The 
HDP quasars from COSMOS, R06 and E94 (red points) do lie at the bottom edge of the 
distribution, but are not as extreme as the two Jiang et al. 'hot- dust-free' quasars. However, 
the L(3.5/zm/L(5100A) ratio alone cannot distinguish the HDP from the galaxy -dominated 
SED shape (e.g., small black points), because they have similar rest frame 3.5 /im to 5100A 
luminosity ratios. 

We plot the 6 quasars from Jiang et al. (2010) on the mixing diagram (Figure [51 magenta 
dots). We find that three of the 6 Jiang et al. (2010) z ~ 6 quasars (J0005-0006, J0303-0019 
and J1411+1217) qualify as HDP quasars. This ~ 50% HDP fraction at z ~ 6 suggests, at 
low significance, that HDP quasars fraction continue to grow with redshift. All three are 
class II, i.e. their NIR emission may be a continuation of the accretion disk. J1411+1217, 
which was classified as a normal quasar in Jiang et al. (2010), is actually HDP by our 
criteria. 

One difficulty with i/L ll (3.5/j,m) /i/L v (51Q0A) is that this ratio does not allow for a range 
in optical slope, such as would be induced by moderate reddening. The slope for the rest 
frame optical SED of J0303-0019 (ccopt ~ 0.6) is flatter than a standard accretion disk, so 
that a continuation of the power law to the rest frame infrared underestimates the hot dust 
contribution. De-reddening the optical SED leads to a requirement for a host dust excess. 
So this quasar should be counted as "hot-dust-poor" rather than "hot-dust-free". Another 
difficulty with z/L I/ (3.5 / um) / uL^^WOA) is that sources with weak hot dust emission but some 
cool dust emission are not selected by this parameter (e.g. XMM-COSMOS HDP quasar 
XID=96, Hao et al. 2010). Defining an index based on only two rest-frame wavelength band 
is not so effective at low redshift. The mixing diagram of Figure [5] is more useful in selecting 
the outliers. 



5.4. Origin of the HDP Quasars 

The origin of HDP quasars is unknown. As for the COSMOS HDP quasars discussed 
in Hao et al. (2010), the E94 and R06 HDP quasars are at z < 4, when the universe is more 
than 2 Gyr old. They have more than 1 Gyr from reionization to form a torus. So they 
cannot be the first generation of quasars. They are also luminous enough, and have high 
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enough accretion rates, to support a dusty torus (Elitzur & Ho 2009). 

Either the hot dust is destroyed (dynamically or by radiation), or the dust is not centered 
on the SMBH, i.e., an off-nuclear AGN (Blecha et al. 2011, 2008, Guedes et al. 2010, 
Volonteri & Madau 2008). When an SMBH recoils (or is kicked-out), it is possible bring 
along the adjacent broad line region, but not the dusty torus which, being further out, is less 
tightly bound to the SMBH (Loeb 2007). The dust is thus not centered on the SMBH, i.e. 
an off-nuclear AGN. A good candidate for a recoiling BH has been found in the COSMOS 
sample (Civano et al. 2010), although it is not included in the XMM-COSMOS HDP sample 
due to the galaxy-dominated NIR-OPT SED shape. 

Volonteri & Madau (2008) estimated the cumulative number of off-nuclear AGNs (offset> 
0.2" ) per square degree versus redshift. Their lower limit on the number of off-nuclear AGNs 
is ~ 2 deg~ 2 at z < 1, ~ 8 deg~ 2 at z < 2, and ~ 11 deg~ 2 at z < 3. Considering the 
COSMOS field is 2 deg 2 , the Volonteri & Madau (2008) prediction is consistent with the 
cumulative number of XMM-COSMOS HDP AGNs (2.5 deg~ 2 at z < 1, 8 deg' 2 at z < 2, 
18.5 deg' 2 at z < 3). A detailed comparison between the theoretical distribution and the 
HDP evolution will be reported in a following paper. We could not make similar estimates 
for the R06 and E94 samples, because the survey areas are not well defined. 

Alternatively, the existence of HDP quasars in samples with different selection methods 
and the continuous distribution of dust covering factor imply that the origin of HDP quasars 
could be related to the AGN structure. Misaligned disks will result from discrete isotropic 
accretion events (Volonteri et al. 2007), which will lead to a wide range of covering factors 
(Lawrence & Elvis 2010). For disks with tilt-only warps (i.e. with no rotation of the line of 
nodes), ~ 14% of the type 1 AGNs will have covering factors less than 20% (Lawrence & 
Elvis 2010). This agrees with the HDP fraction in all three samples. 

5.5. Prospects 

The ongoing UKIDSS and WISE surveys will expand the available HDP samples greatly, 
especially in the z > 1.5 range, where evolution can then be sought with high sensitivity. 

The Hao et al. (2010) mixing diagram is a useful tool for selecting non-standard quasar 
SEDs. 
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Fig. 9.— Position in Lacy et al. (2004) color-space of the XMM-COSMOS (square), R06 
(circle) and E94 (triangle) samples. The red solid points are the HDP quasars from these 
samples. Small sets of symbols are for sources with olopt < 0.2, that is galaxy- or reddening- 
dominated sources. 
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Fig. 10. — Luminosity and redshift dependence of the hot-dust abundance introduced in 
Jiang et al. (2010) of the XMM-COSMOS (square), R06 (circle) and E94(triangle) samples. 
The red solid points are the HDP quasars from these samples. The magenta points are the 
three z ~ 6 quasars in Figure 1 of Jiang et al. (2010) that is HDP. Small sets of symbols are 
for sources with o>opt < 0.2, that is galaxy- or reddening-dominated sources. The magenta 
dashed lines are the linear fit and 3a range of the sample in Jiang et al. (2010). The black 
solid lines are the linear fit and 3a range of the XMM-COSMOS, R06 and E94 samples. 



